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Abstract 

We derive a static potential for a heavy quark-antiquark pair propagating in Minkowski time 
at finite temperature, by defining a suitable gauge-invariant Green's function and computing 
it to first non-trivial order in Hard Thermal Loop resummed perturbation theory. The 
resulting Debye-screened potential could be used in models that attempt to describe the 
"melting" of heavy quarkonium at high temperatures. We show, in particular, that the 
potential develops an imaginary part, implying that thermal effects generate a finite width 
for the quarkonium peak in the dilepton production rate. For quarkonium with a very 
heavy constituent mass M, the width can be ignored for T<g^M/127r, where is the 
strong gauge coupling; for a physical case like bottomonium, it could become important 
at temperatures as low as 250 MeV. Finally, we point out that the physics related to the 
finite width originates from the Landau-damping of low- frequency gauge fields, and could be 
studied non-perturbatively by making use of the classical approximation. 
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1. Introduction 



Heavy quarkonium systems, cc and bb, have turned out to provide extremely useful probes for 
QCD phenomenology [I]. On the theoretical side, the existence of a heavy mass scale makes 
these systems more susceptible to analytic treatments than hadrons made out of light quarks 
only, while on the experimental side, decays of heavy quarkonia lead to relatively clean signals 
which are precisely measured by now. It is therefore not surprising that the modifications 
that the quarkonium systems may undergo at high temperatures, are also among the most 
classic observables considered for heavy ion collision experiments [2]. 

Let us recall explicitly the conceptually clean connection that heavy quarkonium provides 
between thermal field theory and heavy ion phenomenology. Consider the production rate 
r^+^- of dileptons (for example, /U pairs) with a total four-momentum Q = + 
from a hot thermal medium. It can be shown that this rate is given by [3] 

- 'l + T^ l-T^ ^M.^Cf(Q), (1.1) 



d4Q 3(27r)5Q2 y ' Q2 Q 

where we assumed > (2m^)^, e is the electromagnetic coupling, tj^i, = diag(H ), and 

(?< is a certain two-point correlator of the electromagnetic current j^{x) in the Heisenberg 
picture, J^{x) = ... + '^ec{x)^^c{x) — ^eb {x)j^b{x): 

roo 



(5f (Q) = /"dty d3xe*Q-^'^J"(0)J''(x)) . (1.2) 



The expectation value refers to (...) = Z^^Tv [eyip(—H /T){...)], where Z is the partition 
function, H is the QCD Hamiltonian operator, and T is the temperature. 

Now, the heavy quark parts in the electromagnetic current induce a certain structure into 
the dilepton production rate. More precisely, they produce a threshold at around ~ (2M)^, 
where M is the heavy quark mass, and a resonance peak (or peaks) near the threshold, 
with a certain height and width. The height and width could in principle be observed, 
as a function of the temperature T that is reached in the collision. Of course there are 
all kinds of practical limitations to this, related to non-equilibrium features, background 
effects, the energy resolution of the detector, etc, but at least some broad features like a total 
disappearance ("melting") of the quarkonium peak should ultimately be visible. 

These circumstances have lead to a great number of studies of quarkonium physics at 
high temperatures. For instance, various types of potential models have been developed [1], 
with non-perturbative input taken from lattice simulations [5]- It is not quite clear, how- 
ever, to which extent potential models are appropriate at finite temperatures, or which 
non-perturbative potentials should be used as input: at asymptotically large distances, the 
standard Polyakov loop correlator is known to fail to reproduce the expected Debye-screened 
potential [6], while many modified descriptions are afflicted by gauge ambiguities [7]. More 



1 



recently, direct lattice determinations of the quarkonium spectral function have been at- 
tempted [8], given that in principle Euclidean data can be analytically continued to Minkowski 
spacetime [9] (though in practice model assumptions need again to be introduced, given the 
very finite number of points in time direction and the statistical nature of the data that are 
available for lattice studies). Finally, similar observables have been considered for strongly- 
coupled TV = 4 Super- Yang- Mills theory, through the AdS / CFT correspondence [TU] . 

All of the studies mentioned, however, either resort to some degree of modelling, or attempt 
to tackle the complete problem by "brute force" on the lattice. It is our philosophy here to 
rather make some more use of the existence of the heavy mass M and the high temperature 
T that characterise the system, and the property of asymptotic freedom of QCD. Indeed, the 
strict weak-coupling expansion does appear to be applicable (within 10-20% accuracy) to hot 
QCD at temperatures as low as a few hundred MeV, once worked out to a high enough order 
and supplemented possibly by numerically determined non-perturbative coefficients, rather 
than complete functions; evidence for this has been obtained through precision studies with 
a large number of independent observables, such as spatial correlation lengths [HI [12], the 
spatial string tension [13l [H] , quark number susceptibilites [ISj [HI , the 't Hooft loop 
tension [181 [H] ) a^id perhaps also the equation-of-state [20] . 

More precisely, our goal here is to derive, through a resummed perturbative computation, a 
static potential that a heavy quark-antiquark pair propagating in Minkowski time at a finite 
temperature feels. This relatively simple computation illustrates a few interesting phenomena 
that have not been exhaustively addressed before, as far as we know. Eventually, our goal is 
to take the lessons brought by this study to a practical level, by carrying out a more extensive 
numerical investigation of the properties of heavy quarkonium [21] . 

The outline of this note is as follows. In Sec. [2] we set up the observable to be determined. 
The actual computation is discussed in some detail in Sec. [3l We elaborate on the main 
results in Sec. HI and conclude in Sec. \5\ There are two appendices summarising well-known 
formulae from the literature with our conventions. 

2. Basic setting 

Rather than Eq. ()1.2p . we prefer to concentrate in the following on the correlator 



As we recall in Appendix A, the two time orderings are related to each other by C<((5) = 
ex.p{—(3q^)C^{Q); for q^>2M, C^{Q) is of order unity while C<{Q) is exponentially sup- 
pressed. Furthermore, let us for simplicity set the spatial momentum to zero, q = 0, and 
consider the correlator before taking the Fourier transform with respect to time: 
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Note that translational invariance guarantees that Cy{—t) = C<(t). 

Now, the determination of Cy{Q) around the threshold requires a resummation of the 
perturbative series (this is the case even at zero temperature). A way to implement this is 
to deform the correlator suitably so that we obtain a Schrodinger-type equation, which can 
then be solved "non-perturbatively" , implementing an all-orders resummation. To arrive at 
a Schrodinger equation, we define a point-splitting, by introducing a vector r, and consider 
an extended interpolating operator for the electromagnetic current, rather than a local one 
(we also drop the electromagnetic couplings, and denote c, 6 by a generic quark field ip): 

Cy{t,r)^ yd3x(|(t,x+^)7'^VFr[(t,x+|);(t,x-^)]V;(t,x-^) ^(0,0)7^^0)) • (2.3) 

Here VFr[2;i,xo] is a Wilson line from xq to xi, along a straight path in the direction of r, 
inserted in order to keep the interpolating operator gauge- invariant. Once the solution for 
Cy{t,r) is known, one can return back to the original situation r = 0. We stress that only 
the r = limit corresponds to the physical current-current correlator in Eq. (j2.2p that we 
are interested in. Indeed, the way the point-splitting is carried out may not be unique, but 
these ambiguities should disappear once we set r = in the solution. 

Let us inspect C'>(t,r) in the limit that the heavy quark mass, M, becomes very large. 
Then the heavy quarks are essentially non-relativistic, and should also be weakly interacting. 
Ignoring interactions altogether, a straightforward diagrammatic computation shows that the 
correlator in this limit indeed satisfies a Schrodinger-type equation, 



,a-(2M-| + o(-I,) 



(5>(t,r)=0, (24) 



with the initial condition C'>(0,r) = —6NcS^^\r). Here the factor 6 corresponds physically 
to two heavy spin-1/2 degrees of freedom times a sum over the spatial components of the 
current. For C'<(t, r), an analogous computation in the same limit yields 

C<(t,r)=0, (2.5) 



with the same initial condition; this form is consistent with the symmetry C>(— i) = C<(t). 

Once interactions are switched on, we might expect (apart from the renormalization of 
the parameter M and of the initial condition) to generate a potential, V{r), into Eqs. (j2.4p . 
(j2.5p . appearing in a form familiar from non-relativistic quantum mechanics. To be precise, 
we define the potential to be a function inside the round parentheses in Eq. (12. 4p which scales 
as in the large-M expansion. To simplify the situation a bit, we would like to be able 
to extract the potential directly, without needing to bother about gradients like V^/M. One 
possibility would be to consider the modified object 

(7>(t,r) = (|^(t,^)^M^r[(t,^);(t,-^)]^(t,-^)x 

x|(0,-^)7MW-r[(0,-^);(0,^)]^(0,^)) . (2.6) 
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Here the quarks can be considered infinitely heavy, such that they do not move in space; the 
integral over x has consequently been dropped; and the quark fields at t = point-split. 

To simplify the situation even further, one notes that the infinitely heavy quark propagators 
are themselves represented by Wilson lines, apart from the "trivial" phase factor exp{—iMt). 
Thereby we arrive at a Wilson loop. It should be noted, though, that time ordering still 
plays a role, and the precise specification of the object considered is not as simple as in 
Euclidean field theory. In fact, the Minkowskian Wilson loop is most naturally defined just 
as an appropriate analytic continuation of the Euclidean object. 

To summarise, we can imagine at least two ways of defining what we might call the static 
potential in real time. On one hand, one can start by computing the standard Wilson loop 
in Euclidean spacetime, and then carry out the analytic continuation that leads to the time 
ordering in Eq. (j2.6p . This is the computation that will be described below. On the other 
hand, one can consider directly Eq. (12. 6p . just replacing the quarks by (almost) infinitely 
heavy static ones: this amounts to NRQCD [22j with the omission of terms of order 0{1/M). 
In the heavy-quark limit, Cy{t,r) is dominated by the forward-propagating part, with the 
phase factor ~ exp(— 2iMt); the static potential can then be extracted from an equation 
like Eq. (12. 4|) . We have carried out this computation as well, and do indeed obtain the same 
result as from the analytic continuation of the Euclidean Wilson loop; however, given that the 
NRQCD computation is quite a bit more involved than the Wilson loop one, we concentrate 
on the latter in the following. 

As indicated above, the calculation of a static potential necessarily involves a point splitting 
and thus should be regarded as an intermediate step towards the computation of a current- 
current correlator. This feature is obviously shared by potential model approaches. However, 
we believe that our procedure offers several advantages over the latter. Since our calculation 
is performed in a well-defined field theoretic setting, the connection with the physical real- 
time observable of interest remains manifest, while non-perturbative gauge invariance as well 
as the correct perturbative limit of our static potential are ensured. 

3. Details of the computation 

Since we make use both of Minkowskian and Euclidean metrics, let us start by introducing 
some notation to keep them apart. Minkowskian four- momenta are denoted by capital sym- 
bols, Q, with components g^, while Euclidean ones are denoted by Q, with components q^. 
The indices are kept down in the latter case, and our convention is = {qo,qi) = {%, —<?*). 
Spacetime coordinates are denoted by x,x, and here our convention is = {xo,x^). The 
Euclidean scalar product is thereby naturally defined as x ■ Q = x^q^ = xoqo — x*g*, 
the four-volume integral as J- = JqAxq J d^x, and the thermal sum-integral as f^Q = 
T j'd^q/(27r)'^. Wick rotation amounts to xq ix^, % —iq^. All Matsubara fre- 
quencies we will meet are bosonic: i.e. qo = InnT, n € Z. The temperature is often 
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Figure 1: The graphs contributing to the static potential at 0{g^). Arrows indicate heavy quarks or 
Wilson lines, and wiggiy lines stand for gluons. 

expressed as /3 = T~^. 

3.1. Wilson loop with Euclidean time direction 

Let again W[zi; zq] be a Wilson line from point zq to point zi: 

W[zi; zo] = 1 + ig dx^A^{x) + {igf di^ dy^ A^{S:)A^{y) + . . . , (3.1) 

J za J zn J zn 



where A^^ = A'^T"' and are the Hermitean generators of SU(A'c), normalised as Tr \T"'T \ 
6°"^ /2. The Euclidean correlation function considered is then defined as 



Ce{t, r) ^ ^Tr ( 1^[(0, r); (r, r)] W[{t, r); (r, 0)] W[{t, 0); (0, 0)] VF[(0, 0); (0, r)] 



(3.2) 



where we have for convenience shifted the origin by r/2 with respect to Eq. (|2.6|) . The 
prefactor l/N^ has been inserted as a normalization, guaranteeing that C£;(0,r) = 1. 

We can formally expand Ce in a power series in the coupling constant 5^, understanding 
of course that the infrared problems of finite-temperature field theory necessitate the use 
of resummed propagators in order for this procedure to be valid (cf. Appendix B): Ce = 
+ + ... , where the superscript indicates the power of g appearing as a prefactor. 



The leading order result is trivial, = 1. We now turn to the computation of C"^' . The 
graphs entering at this order are shown in Fig. [H 

The computation of the graphs in Fig. [T] is not quite trivial. The problem is that in order 
to avoid ambiguous expressions (of the type "0/0"), one needs to treat the Matsubara zero- 
modes very carefully. In fact, we treat them separately from the non-zero modes. For the 
0(5^)-contribution to the first graph in Fig. [H for instance, we obtain 



.(2) 



VAAAAAAAAA 
* 



,'/ Cf / dxQ I dyo 
/o JO 



g'CE 



d\ 

(2^ 



^(fo(0,q 
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90^0 
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(3.3) 



where we have inserted the gluon propagator from Eq. (|B.3P : defined Cp = {N^ — l)/2Nc; 
and chosen the vector r to point in the x^-direction: r = (0,0, r). It is perhaps appropriate 
to note that if one would rather attempt to keep all the Matsubara modes together, as an 
expression of the type on the latter row in Eq. (j3.3p . and then try to carry out the sum with 
the usual contour trick \23\ [M| , the zero- mode go = is not treated correctly, because of the 
ambiguity of the expression under r ^ r + n/3, n G Z. 

Summing all the graphs together, the parts proportional to ^, as well as the "longitudinal 
parts" cx q^Qi, from the terms multiplied by (cf. Eq. (|B.2[) ). cancel explicitly. Inserting 
finally the expressions of the projection operators from Eqs. (|B.ip . ()B.2p . we obtain 



d^q e 



iq^r 



+ e 



+ 



(2^)3 
1 



q2+n,(0,q)-+^^^^ 



-iqoT 



^2 



+ q2+nij(go,q) 



+ 



1 



1 



;;2 



+ q2+nT(go,q) 



(3.4) 



We have kept the spatial exponents in a form which makes it clear that the divergent-looking 
term is in fact harmless. 

In order to carry out the remaining sum in Eq. (j3.4p . it is useful to write the resummed 
propagators in a spectral representation. For any function A(zgo) which behaves regularly 
enough at infinity, we can define the spectral density p{q^) through 



1 

2i L' 



A(igo ^ g° + iO+) - Aiiqo ^ q° 



which leads to the inverse transform 

A(igo) 

In particular. 



dgO p(gO) 



vr q^ 



dq^ PEiq'' 



(3.5) 



(3.6) 



(3.7) 



go + + ni?(go,q) J-oo g^ 
and similarly with Ily. The properties of the spectral functions pt,Pe are described in 
Eqs. (fB77l) - (fBl3l) of Appendix B. 

Inserting the spectral representations into Eq. (13. 4p , and making use of Eq. (lA.lSp , we can 
carry out the remaining sums: for < r < /3, 



-iqoT 



8T 
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gO/3 



nB(g°) 
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(3.8) 
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(27rnr)2 + (g0)2 



1 + 6' 



=(/3-T)g» 



where nB{q^) = l/[exp(/3g°) — 1]. Noting furthermore that (cf. Eq. ([3 



vr q 



we arrive at 



q2 + ns(0,q) ' 



(27r)3 



q2+ni5(0,q) 
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ml ]PE{q^.^ + 



4 - IPTlg^q) 



(3.9) 



(3.10) 



(3.11) 



q2 (g0)2y ■- - 

It is useful to note that both the second and the third row in Eq. (|3.11|) are odd in g'^, so 
that the product is even, and the integral could also be restricted to positive values of <^ . 



3.2. Equation of motion in Minkowski time 



Eq. (j3.1ip has a form which can directly be analytically continued to Minkowski time, r — > it. 
Thereby we arrive at the object we are interested in, C>(t,r) (cf. Eqs. (IX2|) . (|AJ|)): 



1 , 
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(3.13) 



We are now in a position to identify the heavy-quark potential. We define it as the quantity 
which determines the time behaviour of our correlator, in accordance with the Schrodinger 
equation: 

idtCy{t,r) = Vy{t,r)C>{t,r) . (3.14) 

Note that this V>(t,r) depends, in general, on the time t and on the time-ordering that is 
used for defining the correlator Cy{t,r). 

In perturbation theory, Eq. (I3.14P needs to be fulfilled order by order. This implies that 



idtC^>\t,r) 
idt&^\t,r) 



0, 

vi^\t,r)C^^\t,r) . 



(3.15) 
(3.16) 
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Eqs. ()3.12p . ()3.13p then immediately lead to the identification 

(2), , 2^ f d3q 2 - e^«3r _ g-igar c 1 



(27r)3 



q2+ni,(0,q) 



+ 



oo ^^0 



+ / — nB(g°)g° 
. vr 



For V^\t,r), the analogue of Eq. (I3l^ reads (cf. Eq. (|23 



(3.17) 



iatC<(t,r) = -y<(t,r)C<(t,r) 



(3.18) 



where the symmetry mentioned below Eq. (|2.2p . or an explicit computation with NRQCD, 



yield yi^^(t,r) = V^">{-t,r) 



(2), 



Discussion of the main results 



We have shown in the previous section that the real-time "Wilson loop" C>(t, r), charac- 
terising the propagation of two infinitely heavy quarks separated by a distance r at a finite 
temperature T, satisfies Eq. (j3.14p . where the potential V^{t,r) is given in Eq. (j3.17p to first 
non-trivial order in . Let us now discuss the behaviour of Vy{t,r) in various limits. 

Perhaps the most interesting limit is what happens at large times, t oo. In this limit, 
the exponential functions exjp{±iq^t) average to zero, unless they are multiplied by a singular 
prefactor, ~ then 



-iqOt 



lim ■ 



27ri 5(g°) . 



(4.1) 



Given that n^{q^) « l/(3q^ for \q^\ <^ T, and that Pe{q^,^) has a term linear in q^, cf. 
Eq. (IB.13p . such a (5-function indeed emerges from the part containing p_B((?'^,q). We obtain 



lim Vi^^ (t, r) 



d'^q 
(2^ 



exp(— rn-DT 



I'Km 



D 



q2 + [5 |q|(q^ + ?ti-d)^ 

ig^TCp 



(4.2) 



(4.3) 



where we have also made use of Eq. (jB.lip . and mD is the Debye mass parameter. The 
r-independent constant in Eq. (j4.3p has been evaluated in dimensional regularization, while 
the dimensionless function (j)(x), 



dz z 

(z2 + 1)2 



1 



sin(zx) 



zx 



(4.4) 



8 



is finite and strictly increasing, with the hmiting values (/>(0) = 0, </>(oo) = 10 

We observe, thus, that apart from a standard Debye-screened static potential, the potential 
y>(i,r) also has an imaginary part. The sign is such that Eq. (13.141) leads to exponential 
decay for C>(t,r). The imaginary part obtained for y<(i,r) has the same magnitude but 
opposite sign, so that Eq. (|3.18p leads to exponential decay as well. We will refer to the 
imaginary part as the "decay width". 

In order to get a rough impression of the numerical orders of magnitude, let us solve the 
time- independent Schrodinger equation with the real part of Eq. (j4.2p : 

^ + ^ _ (,2(^^exp(-mDr) 



M \ dr^ r dr J Airr 



ipir) = E ipir) , (4.5) 



where E denotes the binding energy. If M is very large, g'^CpM/SiT ^ mo, Debye screening 
can be ignored, and we get the regular hydrogen atom solution, with an inverse Bohr radius 
1/ro = g'^CpM/8iT, and a binding energy \Eo\ = Mig^Cp/S-Kf. Treating the imaginary part, 
r, as a perturbation, T ~ (g'^TCp/ Att) J dc^r \ip{r)\'^ (j){mj:,r), it is clear that |r| < g'^TCp/^'K, 
because < (j){m-Qr) < 1; in fact, for romo <C 1 as was the assumption, |r| <C g'^TCp / ^tt. 
We thus note that for a very heavy quark mass M, the decay width can equal the binding 
energy only at temperatures T ^ g'^CpM/Wir. 

Lowering towards more realistic quark masses at a fixed temperature, vq increases, and at 
some point the Debye radius l/mo becomes important. Consequently the binding energy 
decreases fast and, simultaneously, the width T increases, because the wave function spreads 
out, whereby (j){mY:,r) obtains non-zero values. 

Unfortunately, it is non-trivial to make these arguments precise: to fix the scale parameter 
appearing in g^ and m^, another order in the perturbative expansion would be needed; 
once the width is substantial, it can no longer be treated as a perturbation; and the full 
time-dependence of V> {t-,''') should be considered in order to determine the width reliably. 
Nevertheless, for illustration, we may fix the scales inside g^ and by making use of 
another context where several orders are available, namely that of dimensionally reduced 
effective theories [25^ [26| [27] . Concretely, employing simple analytic expressions that can be 
extracted from Ref. [28], 

^'" 91n(9.082T/A^) ' " 31n(7.M7r/A^) ' fo^' = A^f = 3 ; (4.6) 



solving numerically the Schrodinger-equation in Eq. (|4.5p : and treating the imaginary part 
still as a perturbation, we obtain Fig. [2] for the bottomonium system. (For charmonium 
the relevant mass M ~ 1.25 GeV is quite a bit smaller, and moves the patterns to smaller 
temperatures, where our methods are no longer reliable.) We stress that this figure is not 
meant to be quantitatively accurate, but it does illustrate something about the qualitative 
behaviour of the system: namely, that the width increases with temperature; the binding 



^The latter part of the integral can be expressed as Meijer's G-function, with certain arguments. 
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400 



T/MeV 

Figure 2: Very rough estimates for the binding energy and decay width of the bottomonium system 
(M ~ 4.25 GeV), as a function of the temperature, based on Eqs. (|4.2p - (|4.6p . The error bands have 
been obtained by varying from 300 MeV (lower edges) to 450 MeV (upper edges). We expect the 
true errors to be much larger than these bands, for reasons discussed in the text. 

energy decreases; and the width exceeds the binding energy already much before the latter 
disappears (i.e. bottomonium "melts"). 

Let us at this point inspect how the results would change in the classical limit^ For this, 
we need to reintroduce h. After the standard rescaling of the fields, h appears in connection 
with the coupling constant, as g'^h, as well as in the extent of the Euclidean time direction, 
as ph. In the classical limit, therefore, nB{q^) — > l/(3h/f. 

We now note that h cancels from the imaginary part of Eq. (j4.3p . where the combination 
g^h/ (3h = g^T appears. This means that the decay width exists also in the classical limit. 
In fact, nothing but the decay width exists in the classical limit! Moreover, time-ordering 
has no meaning in the classical limit, which explains why the imaginary parts have the same 

(2) (2) 

magnitude in V> (t, r) and V<; (t, r). That a decay width / damping rate should be classical, 
need not be surprising: particularly detailed demonstrations of this have been provided for 
scalar field theory [38j. 

Finally we remark that for T — s- 0, i.e. (3 — > oo, we can restrict to q*^ > in Eq. (j3.17p 
because of the symmetry, and then n-B{q^) — > 0, n'Q{q^) exp(/3g'^) 1, whereby the structure 
of Eq. (jS.lTp simplifies. In particular, there is no imaginary part any more for t oo, only 
the standard Coulomb-potential. 

^ This limit is of course singular in many respects, but can still be formally defined [29], and allows to 
address numerically many important non-perturbative problems in weakly-coupled gauge theories [30]-|35), 
in spite of the presence of significant discretization artifacts [361 137) in the classical lattice gauge theory 
simulations [30] that are employed in the thermal context. 
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5. Conclusions and Outlook 

The purpose of this note has been to present a derivation of a "static potential" that a heavy 
quark-antiquark pair propagating in a thermal medium feels. An integral representation of the 
result is shown in Eq. (13.17p . This potential could then be inserted into Eq. (j2.4p . in order to 
estimate a certain real-time Green's function for the vector current, whose Fourier-transform 
determines the quarkonium contribution to the dilepton production rate. A detailed numer- 
ical evaluation is postponed to future work. 

It is the first term in Eq. (j3.17p which represents the standard time-independent Debye- 
screened potential. We note that this potential behaves as — g'^Ci? exp(— mDr)/47rr as a func- 
tion of r. Thus it differs from the potential extracted from the correlator of two Polyakov 
loops in Euclidean spacetime, which behaves as ~ — 51^ exp(— 2mDr)/47rr at phenomenologi- 
cally relevant temperatures [11], and as ~ — gi^"'""' exp(— mGr)/47rr, n > 0, at asymptotically 
high temperatures [6] , where rriQ, is the lightest glueball mass of three-dimensional pure Yang- 
Mills theory [39], with the gauge coupling gl = g^T[l + 0{g)] [HI [40]. 

Motivated by our perturbative analysis, one can however envisage ways of defining a mod- 
ified (real part of the) static potential, which would be gauge- invariant, possess the correct 
perturbative limit, and allow for a direct non-perturbative measurement with lattice simula- 
tions. Inspecting Eq. (j3.1ip . the correct term is seen to sit on the first row. The corresponding 
functional dependence on r differs from that on the second row in a qualitative way: for in- 
stance, it is linear in r at r <C /3, while the second row is quadratic; and it is non-periodic in 
r — > /3 — r, while the second row is periodic. Both of these functional features can in princi- 
ple also be isolated from non-perturbative data for the Euclidean Wilson loop; whether they 
lead to useful practical recipes for determining the real part of the static potential beyond 
perturbation theory, remains however to be tested. 

Apart from the standard potential, we note that Eq. ()3.17p also has another term, the 
second one, with a fairly rich structure. In the limit t — > 00, the second terms amounts 
to a thermal decay width, induced by Landau damping of the low-frequency gauge fields 
that mediate interactions between the two heavy quarks. The thermal width induces a 
certain width also to the quarkonium peak in the dilepton production rate, thus making 
the peak much wider than at zero temperature. Physically, the Landau damping underlying 
this phenomenon originates from an energy transfer from low-frequency gauge fields that 
are responsible for the static interaction, to the "hard" particles (which have momenta of 
the order of the temperature) existing in the thermal plasma; technically, it originates from 
the "cut contribution" to the gluon spectral function at \q^\ < |q|, cf. Eqs. (IB.12p . (IB.13p . 
Restricting to a very rough estimate, we have shown that this part leads to physically sensible 
qualitative structures, namely a width which increases with the temperature, exceeding the 
binding energy already quite a bit before the latter disappears and quarkonium "melts" . 

It is worth stressing, however, that apart from the limiting value at i — > 00, the potential 
has a non-trivial time dependence, which also plays a role in the solution of Eq. (|2.4p and, 
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consequently, for the quarkonium contribution to the dilepton production rate. 

In order to define the significance of these findings for heavy quarkonium, its finite mass 
M has to be taken into account systematically. For this one can make use of NRQCD |22j . 
allowing to account analytically for all exponential effects like exp(— 2/?M), so that one can 
concentrate on the softer dynamics around the threshold. Perturbatively, this amounts to the 
solution of Eq. (j2.4|) with V> {t, r) inserted inside the round brackets. However, one might 
also be able to go beyond perturbation theory, by making use of the classical approximation 
for the gauge fields appearing in NRQCD, provided that genuinely quantum effects, such as 
the first term in Eq. (j3.17p . are properly taken into account. Hopefully this recipe could be 
formalised such that one would be able to account for the orders {g^fif and {g^hY in the 
"zero-temperature" expansion, but in each case for all orders in the possibly larger finite- 
temperature expansion parameters, g^T/m-£, and g^T/mQ. 
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Appendix A. Real-time Green's functions 

For completeness and to fix the notation, we list in this Appendix some common definitions 
and relations that apply to two-point correlation functions built out of bosonic operators; for 
more details see, e.g., Refs. [23| [24]. 

We use the notation introduced at the beginning of Sec. El with t = and t = xq. 
Arguments of operators denote implicitly whether we are in Minkowskian or Euclidean space- 
time. In particular, Heisenberg-operators are defined as 

0(t, x) = e'"^6{0, x)e-'"^ , 0(r, x) = e^^O(0, x)e-^^ . (A.l) 

The thermal ensemble is defined by the density matrix p = ex.p{—(3H). We denote the 
operators which appear in the two-point functions by (j)^{x), (jf{x)] they could either be 
elementary field operators or composite operators. 

We can now define various classes of correlation functions. The "physical" correlators are 
defined as 



C^^iQ) = j dtd^xe'^■''(^4>^'ix)^''iO)^ , (A.2) 
Cf(Q) = Jdtd^^e'^-''(^4>''{0)4>^'{x)'j , (A.3) 
pM-(Q) ^ fdtd^^e'<^--(^^[r{x),r{0)]) , (A.4) 
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iQ-x 



(A.5) 



where p'^" is called the spectral function, while the "retarded" / "advanced" correlators can 
be defined as 



C'j^^(Q) = ijdtd 
C^^iQ) = ijdtd 



r{x),rio)\9it)) , 



(A.6) 
(A.7) 



On the other hand, from the computational point of view one is often faced with "time- 
ordered" correlation functions, 

^^(Q) = jdtd^yie'^-''(^^{x)4>''{0)e{t) + 4>''{0)4>^{x)e{-t)y (A.8) 

which appear in time-dependent perturbation theory, or with the "Euclidean" correlator 

C^'iQ) ^ jyTjd'^e"^-^{rixWiO)), (A.9) 

which appears in non-perturbative formulations. Note that the Euclidean correlator is also 
time-ordered by definition, and can be computed with Euclidean functional integrals. 

Now, all of the correlation functions defined can be related to each other. In particular, 
all correlators can be expressed in terms of the spectral function, which in turn can be 
determined as a certain analytic continuation of the Euclidean correlator. In order to do 
this, we may first insert sets of energy eigenstates, to obtain the Fourier-space version of the 
so-called Kubo-Martin-Schwinger (KMS) relation: C'^''{Q) = e-^«°C'^''(Q). Then pi^^iQ) = 
[C^'^iQ) - C'^{Q)]/2 and, conversely. 



(Q) = 2[1 + n^{q'')]r{Q) , &<{Q) = 2nB{q')nQ) ■ 
Moreover, Af^^{Q) = [1 + 2nB{q^)]p^^ (Q)- Inserting the representation 



r 

e{t) = i 

J—r 



,2Tr L0 + iO+ 



(A.10) 



(A.ll) 



into the definitions of Cr, Ca, we obtain 

do; p^''{ui,q) 



TT LO — — zO+ 



/oo J 
-oo ^ 



duj p^"'{u},q) 



Doing the same with Ct and making use of 

1 



A±iO+ 



p(i)Tm<5(A), 



(A.12) 



(A.13) 
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produces 



+ 2p'^-(gO,q)nB(g"). 



(A.14) 



J vr q'^ — uj + iO+ 

Finally, writing the argument inside the r-integration in Eq. ()A.9p as a Wick rotation of the 
integrand in Eq. (|A.2p . which in turn is expressed as an inverse Fourier transform of {Q)^ 
for which Eq. (|A.10p is inserted, and changing orders of integration, we get 



du) p'^'^{uj, q) 



'0 J-oo ^" J-oo ^ ^-iQo 

This relation can formally be inverted by making use of Eq. (lA.lSh . 

P^"('?°,q) = ^[C'r(-ib° + ^0+],q) - C>^''i-z[q' - iO+],q)" 
We also recall that bosonic Matsubara sums can be carried out through 



90 



+ ^2 



{cdr + d)TY, 



go 



2E 



cE + d)e^^-^^^ + {cE + d)e^' 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



where qq = lirnT, with n an integer, and we assumed < r < /?. This equation can be used 
as a starting point for determining the sums in Eqs. p.Sp . (j3.9p . 



Appendix B. Resummed gluon propagator 

Introducing the projection operators (see, e.g., Refs. 



P^iQ) 



the Euclidean gluon propagator can be written as 



Q 



PLiQ) 



.Q2 + nr(Q) Q^ + Ue{Q) (Q2)2 
where ^ is the gauge parameter. The Hard Thermal Loop \^T\ H2] contributions read 



(B.l) 
(B.2) 

(B.3) 



nT(0) 



2 \ q2 2|q| 



1 



im) 

q2 



In 



m + |q| 



mo 



1 



im) 

q2 



1 



^90 - |q| 
m m + |q| 

2|q| iqo - |q| 



(B.4) 
(B.5) 
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where % denotes bosonic Matsubara frequencies, and 



mf, = q 1 

° ^ 13 6/ 



(B.6) 



In the hmit 1% — s- but with |q| 7^ 0, Ht 0, II^; — > m^, while for |q| with Ojq ^ 0, 
Hr, He ^ my3. 

After analytic continuation, 1% ^ + iO~^ , the propagators become 

1 1 



Q2 + nr(i,)(go, q) -('?° + ^0+)2 + q2 + nr(i5)(-i(<70 + iO+),q) 



where 



nTH(g° + io+),q) 

n£;(-i(g° + iO+),q) 



f (q^)"^ q^ 



q- 2|q| 



0\2 



In 



g° + iO+ + |q| 
g° + iO+ - |q| 



0\2 



f , q° + iO+ + Iql 
■ m 



2|q| gO + ^0^ 



(B.7) 

(B.8) 
(B.9) 



For \q^\ > |q|, IIt,IIe are real. For \q^\ < |q|, they have an imaginary part. In particular, 
for \q^\ <^ |q|, we get 



2|q| 



+ 2 



,0\2 



+ 



n, 



mj3< 1 + in 



q- 



2|q| 



+ ... 



(B.IO) 
(B.ll) 



The Hard Thermal Loop resummation of course only describes the behaviour correctly for 
|q| ~ gT; for instance, for |q| ~ g^T, loop corrections within the Hard Thermal Loop effective 
theory are large, and the small-frequency behaviour of the full H^ is determined by a "colour 
conductivity" rather than Eq. (|RT0]) : Ht ~ -icrq°, where a ~ r/ln(l/5) [351 [ii]. 

The spectral functions pt,Pe that are needed in the text follow by taking the imaginary 
part, or discontinuity (cf. Eq. (j3.5p ). of the right-hand side of Eq. (|B.7p . In particular, 



Pr(9°,q) 



/5i?(g°,q) 



TTsign{q^)5\^{q^y - q - ReHr 
2 9° 



' 7rsign(g0)5(^(g0)2 -q2 -ReHe) 



2|q|(q2 + m^) 



2 \2 ' 



\q \ > |q| 



\q \ <■ |ql 



\q \ > |q| 



Irl < |ql 



(B.12) 



(B.13) 



Again, "soft" loop corrections to these expressions are large for ultrasoft momenta, Iq[ ~ g^T. 
Moreover, loop corrections are also significant around the plasmon poles, i.e. the (5-functions 
in Eqs. (IB.12p . (jB.lSh . where a finite plasmon decay width gets generated f35] . 
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